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Abstract 

Following Blumenthal et al. [1] and Piran et al. [2] I discuss a simple model for 
the gravitational formation of voids. This model enables us to estimate the primordial 
perturbation required on the scale of voids. A comparison of this condition with the 
constraints on the spectrum of fluctuations derived from CMBR observations reveals that 
gravitational growth of large voids require more power on the scale of ~ 1500km/s than is 
generally accepted in CDM models. For a generic power law primordial spectrum P{k) oc 
A;” gravitational growth of voids requires n ^ 1.25, which is marginally compatible with 
the observed CMBR fluctuations. If future measurements will show that the amplitude 
of fluctuations required to produce the voids is too large they will imply that the voids 
did not grow gravitationally and that galaxies do not trace matter on very large scale. 

To appear in the Proceedings of the XXXth RENCONTRES DE MORIOND Workshop on 

Dark Matter in Cosmology 


1 INTRODUCTION 

Recent observations of the galaxy distribution strongly suggest that galaxies tend to lie on 
wall-like features which bound large empty regions, forming a closely packed network of voids 
[3, 4, 5, 6, 7]. This property seems to be generic since it is present in all the slices in the 
northern and southern hemisphere that have been completed by the ongoing CfA2 and SSRS2 
surveys [6, 8]. Further evidence for this picture has been provided by the red-shift maps of 
the deep survey being conducted by [9]. Combined these surveys are qualitatively suggestive 



of a void-filled Universe. Gross estimates of the typical size and under-density suggest that the 
voids have a characteristic scale ranging from 2500 to 5000 km/s in diameter with a typical 
under-density estimated at about 20% of the mean (see [If], for the Bootes void). Currently 
there is no clear statistical measure for the properties of voids. Fig. 1 depicts a preliminary 
distribution of voids in the SSRS2 survey obtained using a new void search algorithm by El-Ad 
et al. [10]. The hgure supports the qualitative picture described above. 



Fig. 1: Voids in the SSRS2 survey [8]. The survey is in the south galactic hemisphere and its 
boundaries are: 20.Oh < RA < A.GGhours and declination —40° < 8 < —2.5°. The figure 
depicts a .slice of constant declination in the range —12.5° < 8 < —7° from [10] 

Following Blumenthal et al. [1] Dubinski et al. [12] and Piran et al. [2] I describe here 
a simple model for the gravitational formation of voids. According to this model voids grow 
gravitationally from negative density fluctuations which are now at the stage of shell crossing. 
The condition of shell crossing constraint the primordial fluctuation spectrum on the scale 
of voids. This criterion enables us to incorporate the existence of voids in the quest for the 
primordial spectrum of fluctuations. Alternatively, if one can show that this constraints cannot 
be satished by the primordial fluctuation spectrum we will have to conclude that light does not 
trace matter on the scale of voids. 





















































































2 GRAVITATIONAL GROWTH OF VOIDS 


Consider, first, the growth of a single isolated spherical inverted top-hat void in an 0 = 1 
universe. We assume that at some initial time C' (e-g- at horizon crossing) there is an inverted 
top-hat density distribution; i.e., for r < i?p p = pc(l ~ q) where pc is the critical density and 
Ci = 1 — Op and for r > Ri, p = p^- Initially this perturbation grows linearly: the under-density 
increases with time while its size grows only with the expansion of the universe. Eventually the 
perturbation becomes non-linear and its size increases faster than the expansion of the universe. 
Finally at a time given by €,{tsc/t^y^^ sa Scrzt ~ 2.5, (i.e. when £*(1 -g Zi)l{l + z,c) ~ Scr^t) 
shell crossing begins (with the shell just inside the boundary being the hrst to cross). We dehne 
the instant of shell crossing as the moment at which the void forms. This dehnition might seem 
somewhat arbitrary. However, at Cc the co-moving radius has grown and it is 1.7 times 
the initial co-moving radius Ri. The density in the void is then 1/1.7^ ~ 0.2 of the average 
density, which is the frequently quoted value for the galaxy under-density in the observed voids. 
The requirement for shell crossing may perhaps be too extreme, although it naturally leads to 
under-density values which are roughly comparable to observations. Relaxing this condition 
and requiring the under-density to be 30% of the mean (which corresponds to an expansion over 
the co-moving scale by a factor of 1.45) implies that the amplitude of the initial perturbation 
could be decreased by a factor of 1.5. 

Despite the simplicity of the adopted model, its main features are supported by N-body 
experiments of simple conhgurations of several interacting voids carried out by Dubinski et al. 
[12] . In these simulations it was found that the condition for formation of an isolated void 
holds approximately, even for highly non-spherical perturbations and when other negative or 
positive perturbations are present. These numerical simulations also reveal that voids at shell¬ 
crossing are the most prominent being delineated by high-density contrast walls. Only tenuous 
traces of smaller scale voids, past the shell-crossing phase, are seen producing a distribution 
that roughly resembles that actually observed for the galaxies. 

Two other properties are critical to our model. First, the largest voids at any epoch are 
those reaching shell crossing at that epoch, since after shell-crossing voids grow relatively slowly. 
Second, if the fluctuations are Gaussian-distributed then when rms negative fluctuations on a 
given scale reach shell crossing the corresponding voids will also occupy a large fraction of the 
volume of the universe, resembling a void-hlled universe [12]. We use this fact to identify e with 
(7, the rms mass fluctuation on the scale of voids. It might be that in a more realistic model 
1.5(7 or 2(7 fluctuations will be sufficient to produce a void-hlled universe. This introduces an 
additional uncertainty in our model which is difficult to estimate. 

The quantity Ciil + Zi) jil + z) equals to the amplitude of the perturbation 6M/Ml had it 
continued to grow linearly until today. It can be expressed in terms of the linear perturbation 
spectrum Phik) as: 

{(¥)/«') • <'' 

yV{kR) is a window function, for which we use the Fourier transform of a Gaussian which has 
been normalized to have the same spatial volume as a top hat of radius R\ 

Wix) = exp[—-a;^(2/97r)^/^] . (2) 

This choice of a window function assures the convergence of the integrals in Eq. 1. For a given 
linear power spectrum we can equate 8MjM to Scrit ■ This provides an equation for R the scale 
of voids today. Alternatively, if the scale of voids is known then then this equality becomes an 
integral equation for PL[k). 



3 LIMITS ON THE POWER-SPECTRUM FROM CMBR DATA 


The same mass fluctuations that determine the gravitational growth of voids are also respon¬ 
sible for the CMBR fluctuations measured today. For a given form of the power-spectrum we 
normalize the power-spectrum using the observed CMBR fluctuations. This will determine the 
co-moving scale that satisfies the condition SMjM = Scrit and in turn will predict the size of 
voids today, consistent with the CMBR fluctuations. This prediction should be compared with 
the observed scale of the voids. 

The fluctuations in the CMBR on scales larger than one degree, which are of interest here, 
are dominated by the Sachs-Wolfe effect [13]. In a flat (0 = 1) universe, C{6) can be expressed 
in terms of the current linear power spectrum, PL[k), as: 


C{9) _ sinfc( 
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where 

( = 2Rhsin{0/2) , (4) 

and Rfi = 2clH q is the horizon size. 

For a given C'(0), Fq. 3 is another integral equation on PL{k). For a given functional shape 
for PL{k)^ such as the case for a CDM model, Fq. 3 determines the amplitude of the linear 
power spectrum PL{k). 


4 VOIDS AND CMBR ANISOTROPY 


Using the results of the preceding sections we can now examine the size of typical voids today 
which are consistent with the constrains on the power-spectrum imposed by the CMBR experi¬ 
ments. Following Blumenthal et al. [1] and Piran et al. [2] we choose a power spectrum, 
of the form PL{k) = Ank^. This might be justified since we are interested in extrapolation from 
the scale of COBF measurements to the scale of the voids. We discuss, later, qualitatively the 
implication of replacing this form by the linear CDM power spectrum. For each n we integrate 
Fq. 3 and define, Fn{6) such that : 
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Using the COBF data [14] we then determine the appropriate normalization constant. An and 
integrate Fq. 1 to obtain an estimate of 6M/Ml as a function of the power index n: 
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We have used the velocity, v = FIqR^ instead of the radius, i?, to obtain a formula that is 
independent of the Hubble constant. We also remind the reader that since at shell crossing 
Rsc ~ 1.7Ri^ a typical void with a diameter of 2Rsc today corresponds to an initial perturbation 
with a co-moving radius of Ri. 

The condition imposed by the existence of voids that evolve gravitationally from rms fluctu¬ 
ations can be expressed by equating Fq. 6 to the value of Scrit = 2-7 discussed in section II {Scrit 
could be as low as 1.8 (if we consider under-density of 30% in the voids or if 1.5a rather than 
la fluctuations produce the typical voids). Piran et al. [2] calculated SM/Ml as a function of 



Vf = 3.4 X u, where v corresponds to the co-moving size in km/s of the unevolved primordial 
perturbation for different values of n. These results are summarized in Fig. 2 and in Table 
1. They hud that in order to form typical voids of 5000 km/s in diameter we must consider a 
power-spectrum with n sa 1.25. This could be the case if the primordial spectrum is steeper 
than expected from inflation. 
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Fig. 2: Predicted 6 M/Ml vs. the void diameter (measured in km/s ) for various power laws 

spectrum, whose amplitude fits COBE, from [2]. 

Table 1. I Comparison of CMBR and voids for different spectral indices [2]. 
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1) 6M/Ml on a scale that corresponds to voids with a diameter of n = 5000A;m/sec today. 

2) Diameter of the voids (using 6M/Ml = 2.7) in km/sec. 

3) \JC (0) convolved with a FWHM beam of 1.2°. This should be compared with the 2(7 
upper limit of Gaier et al. [15] of 1.4 • 10“®. 

For a Harrison-Zel’dovich spectrum [n = 1), consistent with the COBF observations and 
the South Pole anisotropy upper-limit, 6M/Ml is sa 1.2 at the scale of 5000 km/s , and thus 
such voids are not expected to be typical. The typical diameter of volume-hlling voids should 
be in this case sa 3500 km/s . Voids as large as 5000 km/s could still be formed from 2(7 
fluctuations and they would represent the high-end tail of the distribution of void sizes. 

Of course, some of the difficulties of forming large voids could be alleviated if we consider 
the possibility that voids are essentially formed when the under-density is of the order of 30%. 
In this case the required amplitudes are decreased by a factor of 1.5 and voids with diameters of 






4200 km/s can form from f cr fluctuations for Scrit =2.7, (3500 km/s for Scrit =4.5) for a Harrison- 
Zel’dovich spectrum. We note, however, that in this case we would not expect them to hll in 
the volume. We also expect the voids to have shallower prohles and less sharp boundaries. 
Unfortunately, the available data is far from adequate to address these details. 

Alternatively, if 1.5a fluctuations are sufficient to produce a volume hlled Universe, (or if 
the voids are not volume hlling) this will also reduce by a similar factor the required amplitude 
for formation of voids and 4200 km/s voids could form with a n = 1 spectrum. However, one 
has to relax both assumptions - i.e. to require that the observed voids are shallower and that 
they are produced by 1.5a events - to explain the gravitational formation of 5000 km/s voids 
with such a spectrum. 

Piran et al. [2] have considered only generic power-law spectra, but this approach can also 
be used to predict the size of typical voids for any specihc model of primordial fluctuations. If 
we consider CDM models we should replace the power-law spectrum by another function, such 
as the one proposed by Efstathiou, et al. [17]: 
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where 0 is the closure parameter and h the Hubble constant measured in lHHkm j.secjMpc. This 
function has the same power on large scale as a pure power law but less power on smaller scales. 
Since CODE measurements are on large scale, substitution of Eq. 7 to Eq. 3 will give effectively 
the same normalization. However, its substitution to Eq. 1, will give lower 6M/Ml on a given 
scale R. Hence the overall effect of replacing the generic power law spectrum by a more realistic 
CDM model is to increase the required power law (which appears also in the CDM perturbation 
spectrum. Recent calculations along this line by Lecar and Gregory [16] show that the required 
index increases to n = 1.5 (with the standard assumptions described above). This result is 
mildly in conflict with smaller scale (sa U) CMBR measurements, whose comparison with the 
CODE data suggest that 1 < n < 1.25 [18]. With the standard parameters of CDM: a Harrison- 
Zel’dovich spectrum n = 1, h = 0.65, and the CDM transfer function (Eq. 7), normalized by 
the two-year CODE results, they shows that the amplitude required to form 5000km/s voids 
gravitationally is achieved only for 3-sigma events. 

Lecar & Gregory [16] have also examined an open Universe with cosmological constant such 
that 0 = 1 but flmatter < 1- In this universe if/^, (which is proportional to 0“^) increases 
and this results in a much larger separation between the scale of CODE (> 10°) and the scale 
of voids (5000A;m/sec). At the same time the perturbations gradually stop growing (after 
1 A z ^ The combined result of both effects is that the required power law index, n, 

remains practically unchanged. The examination of a 0 < 1 universe is much more complicated 
(see for discussion of CMBR fluctuations [19]) and it is under consideration now. 


5 DISCUSSION 

I have presented here a quantitative model for the incorporation of the information on the 
existence of voids to the analysis of the primordial perturbation spectrum. Our basic result is 
that the “common” existence of voids on scales ranging from 2500 to 5000 km/s requires, if 
they formed gravitationally, more power in the range 800 to 1500 km/s than predicted from 
the standard unbiased CDM model. This implies that non-linear effects should be important in 
the evolution of perturbations on these scales. This contrasts with the conventional idea that 
scales beyond about 8 h~^ Mpc are still in the linear regime. With our standard assumptions: 
Scrit =2.7, Icr fluctuations produce the observed voids ( with an under-density of 20%) and 



P{k) oc A;”, a power law index of n = 1.25 is required to produced a 5000 km/s void filled 
Universe and that a Harrison-Zel’dovich, [n = 1), spectrum can produce gravitationally, only 
3500 km/s voids. Similar considerations constrain the largest possible void to about 6000 
km/s in diameter. Voids as large as those suggested by Broadhurst et al. [20] cannot form 
gravitationally and they cannot correspond to voids in the matter distribution. 

We expect that in our model there will be a relationship between the abundance of voids, 
whose origin is negative density fluctuations and rich clusters that originate from positive ones. 
For a void-hlled universe with a scale of 5000 km/s we expect (assuming that rich cluster form 
from positive density fluctuations with the same initial amplitude as the negative perturbations 
that form voids) the abundance of dark matter clusters with masses of 5 X Mq to be of 
the order of 8 X 10“®i7“^/Mpc^, corresponding roughly to galaxy clusters of 5 X lO^‘^Af 0 . This 
should be compared with an Abell richness 1 cluster density of one per Mpc)^ and with 

the density of Abell clusters of richness 2, one per (95/i“^ Afpc)^ ([21])- It seems that our model 
predicts more rich clusters than observed. However, given all the uncertainties involved in 
the model and the data available for clusters we cannot discard the possibility of non-linearity 
on the basis of the existing data on cluster abundance. Specihcally, if voids form from 1.5a 
fluctuations the symmetry between positive and negative fluctuations is broken, reducing the 
predicted cluster density. Alternatively, it is possible that the primordial perturbation spectrum 
is not Gaussian. 

The basic assumption of our model is that the “voidy” nature of the galaxy distribution 
is a natural consequence of negative rms amplitude fluctuations which grow in size during 
the evolution of the universe producing an essentially void-hlled distribution. Although the 
details of our model are uncertain and the information on voids is still sketchy, our formalism 
can be used to quantitatively assess the likelihood of a given model of structure formation to 
generate voids. The underlying assumption of our approach is that voids grow gravitationally, 
that light traces matter at large scales and therefore the large voids in the galaxy distribution 
correspond to real voids in the matter distribution. If further evidence will suggest that there is 
indeed contradiction between the power required for gravitational formation of voids and other 
observations this will imply that the voids in the galaxy distribution do not trace the matter 
distribution. This will further imply that the observed power spectrum measured from large 
scale galaxy surveys does not reflect the underlying matter power spectrum on the scale of the 
voids. 

It is important to stress that the simple model in which galaxies do not form in under- 
dense regions, and consequently, the under-density in the the dark matter is greatly enhanced 
in the galaxy distribution is not sufficient to resolve the large voids problem. In this model 
we would expect that within the observed voids there will appear marked ’’walls” of galaxies, 
which reflect the over-density on smaller scales that exists in the dark matter. Put differently, 
in order for a biasing scheme of this kind to create a large void which does not follow the dark 
matter distribution, it will also have to erase structures corresponding to large amplitude dark 
matter perturbations within the void. A simple interpretation of this result is that the voids 
do not result from a gravitational phenomenon. 

In summary there are two potential implications to our analysis. If this simple model is valid 
then either the power spectrum of primordial perturbation has a slope n ^ 1.25 or the voids in 
the galaxy distribution do not result from a gravitational phenomenon. The former possibility 
is incompatible with simple inflationary models (but it is possible with more complicated ones 
that allow for two scalar helds). The second possibility means that galaxy distribution on scales 
of ^ 3000km/s does not follow the dark matter distribution. 
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